Mott insulators and correlated superfluids in ultracold Bose-Fermi mixtures 
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We study the effects of interaction between bosons and fermions in a Bose-Fermi mixtures loaded 
in an optical lattice. We concentrate on the destruction of a bosonic Mott phase driven by repulsive 
interaction between bosons and fermions. Once the Mott phase is destroyed, the system enters a 
superfluid phase where the movements of bosons and fermions are correlated. We show that this 
phase has simultaneously correlations reminiscent of a conventional superfluid and of a pseudo-spin 
density wave order. 
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Cold atoms loaded on optical lattices provide a new 
tool to study condensed matter systems. These systems 
are experimental realizations of lattice models for quan- 
tum interacting particles such as Hubbard models. They 
provide a simple and controllable way to study various 
interesting phenomena such as the superfluid to Mott in- 
sulator (MI) transition [l[ or Anderson localization 
Recently, there has been growing interest in the physics 
of boson-fermion mixtures at low temperature. These 
mixtures served originally to lower the temperature of 
the fermions by sympathetic cooling Q with the bosons 
in order to study the low temperature physics of inter- 
acting fermions, such as superfluidity in the BCS or BEC 
regime Q. 

In some cases the interactions between fermions and 
bosons cannot be neglected. Fermions have, for exam- 
ple, been shown to reduce the phase coherence of weakly 
interacting superfluid bosons Theoretical studies 

d, 0, [1] have suggested the possibility for the system to 
form pairs (polarons) of bosons and fermions (particle- 
particle or particle- hole), with these composite fermionic 
pairs then adopting different phases such as normal Fermi 
liquid or solid phases. This system was studied numer- 
ically in one dimension by Sengupta et Pryadko Q and 
shows a commensurate filling phase where the sum of 
the number of bosons, TVb, and fermions, iVf, is equal to 
the number of lattice sites L: TVb + Nf = L. This phase 
has been studied simultaneously in the special case where 
Ny, = iVf = L/2 by Pollet et al. who showed that 
it could be understood, in the strong interaction limit, 
as the pseudo-spin density wave (PSDW) phase previ- 
ously proposed by Kuklov and Svistunov [6|]. This kind 
of system has also been studied in the presence of har- 
monic traps [m , but the comparison with experiment 0| 
is problematic. 

In this letter we study the effect of added fermions on a 
bosonic MI phase [H at low temperature. As the boson- 
fermion interaction is increased, the MI is destroyed and 
the system enters a superfluid phase (SF). We will re- 
view the properties of this phase and show that it can be 
crudely described as composed of superfluid bosons and 
of the aboved mentionned PSDW phase. In addition, we 




FIG. 1: At low Ubf, the system consists of free fermions on a 
bosonic Mott insulator substrate. As Um becomes larger the 
Mott insulator is destroyed and the system becomes superfluid 
(SF). There is a crossover between two different superfluid 
regions (see the text). Finally for large values of Uh^, fermions 
and bosons phase separate (PS). 



study the properties of the PSDW phase when there is 
an imbalance in the boson and fermion populations but 
with Nb + N{ = L. 

The system is described by a one-dimensional repulsive 
Hubbard model 
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where the operator b] {bi) creates (destroys) a boson on 
site i while fj and fi are the corresponding fermionic 
operators. The first term in ([T]) describes the hopping 
of bosons and fermions between neighboring sites. We 
have chosen the hopping amplitude t = 1 , which sets the 
energy scale, to be the same for both species, (rt-) is 
the bosonic (fermionic) number operator at site i and Uhh 
and Uhf are the boson-boson and boson-fermion contact 
interaction terms. 

The bosonic, G-t,{x) = (6-6^+:^), fermionic, G[{x) = 
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FIG. 2: Schematic representation of the destruction of the 
bosonic Mott insulator. As Um is increased, the sites that 
were occupied by a boson and a fermion are replaced by sites 
occupied by two bosons. The number of bosons per site is no 
longer 1 and the system can become superfluid. 
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ififi+x), and composite, Gbf(a;) = {bj fl^^bi+^, f^) , Green 
functions, give access to the phase correlations of bosons, 
fermions and particle-hole pairs, respectively. A non 
vanishing bosonic superfluid density pb = (W^b}/2/3 or 
a fcrmionic stiffness pf = (Wf^)/2/3 shows the presence 
of slowly (algebraically) decaying phase correlations {Wh 
and Wf being the bosonic and fermionic winding num- 
bers). The paired (-I-) and counter- rotating (-) super- 
fluid densities, p^f = ((Wb ± W{)'^)/2f3 are both equal to 
Ph + Pf if the movements of the two types of particles are 
uncorrelated and differ if the movements are correlated 
(or anticorrelated). 

We study this model using a canonical quantum Monte 
Carlo method based on "worm" updates which is ex- 
tremely efficient for the measurements of equal time 
Green functions ■ This method was modified to tackle 
the case of mixtures by allowing simultaneous moves of 
the bosons and fermions [ist . 

We set the number of bosons A^b to be equal to the 
number of sites L = 64 and set Uhh = 8, to obtain a solid 
Mott phase for the bosons in absence of interactions with 
the fermions. We studied the evolution of the system as 
Uhi is gradually increased for two different values of the 
number of fermions (A^f = 15 and 35). 

For Uhf < 6.5 the superfluid densities (Fig. [1]) show 
that the system behaves as in the non interacting limit 
(MI-I-FF). The bosons form a MI where ph = 0, each site 
is occupied by only one boson and jumps are forbidden 
because they cost a large energy Uhh- On the contrary, 
the fermions move freely as they experience the same in- 
teraction energy Uhf on all the sites {pf ^ 0). When Um 
becomes larger than Uhh the Mott insulator is destroyed 
and the bosons become superfluid. In this case, it is en- 
ergetically favorable to form pairs of bosons on the same 
site and to have isolated fermions (Fig. [2]). The system 
is then split into doped superfluid bosonic regions and 
fermions. If the repulsion between bosons and fermions 
is not too strong (6.5 < Uhi < 10), the bosons can pass 
from one superfluid region to the next one by the tun- 
nelling through fermion-occupied sites and a superfluid 



FIG. 3: Finite size scaling of different densities for two values 
of the interaction. The SFa regime observed in Fig. [T] around 
Uhf — 7 disappears for larger size whereas the SFb regime 
persists for Uhf — 8. 

(SF) phase is established. For Uhi > 10 the tunnelling is 
suppressed and there is a phase separation between in- 
dependent superfluid regions separated by fermionic re- 
gions. This case is difficult to study numerically because 
the system is frozen. 

When the MI is destroyed, two distinct regions [l^ . 
separated by a crossover, appear: a first region (SFa, 
6.5 < Uhf < 7.5) where correlations of the movements 
of bosons and fermions are not visible {ph + Pf = Phf) 
and a second one (SFb, 7.5 < Uhf < 10) where p^^ > 
pj^f . In this latter region moving a fermion independently 
of the surrounding bosons costs Uhf, while exchanging a 
boson and a fermion has no energy gap associated. These 
favored exchanges lead to the observed anticorrelation of 
winding numbers. 

The SFa region disappears in the thermodynamic limit 
(or at least becomes very narrow) as can be observed in 
Fig.[3l where ph goes to zero with 1/L for Uhf = 7. On the 
contrary, around Uhf = 8, the SFb phase persists in the 
large size limit. We then expect a direct transition from 
the MI-I-FF phase to the SFb regime (Fig. [3] corresponds 
to a density of fermions Nf/L ~ 0.55 but the two values 
of Uhf correspond to the same phases as in Fig. [1] where 
Nf/L ~ 0.23). 

In Fig. m we show the Green functions in the MM-I-FF 
(left) and in the SFb (right) phases. In the Alott 
phase, Gh{x) decays exponentially while Gf{x) decays 
algebraically following the free fermion Green function 
Gff{x) = sm{NfTTx/ L)/ sm{7Tx/ L). The two species are 
independent, in spite of the large value of Uhf, and con- 
sequently Ghf = GhG[. 

In the SFb, all the Green functions decay algebraically. 
The slowest decay is obtained for the bosons with Gh oc 
sin(7rx/L)~°-^^ while Gf oc sin(Aff7ra;/L) • sin(7rx/L)~^-^. 
In this case, Ghf is much larger than both the product 
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FIG. 4: Left : For C/bf = 5 the bosons are in a Mott insulating 
state and Gb decays exponentially. The fermions follow the 
free fermions behavior. Right: In SFb regime, we observe 
phase coherence in all the cases. Gm is larger than G{. 



GhG{ and Gi and follows (approximately) Gff. This be- 
havior confirms that boson-fcrmion pairs arc important 
degrees of freedom in this system (in the SFa region, on 
the contrary, G{ is larger than Gbf fl^V 

Returning to Fig. [21 an intuitive but crude physical in- 
terpretation is, that our system can be roughly described 
as composed of A'f supcrfluid bosons and a layer com- 
posed of N{ fermions and L — Nf bosons. For such a 
"filled" layer, the description of the system as pseudo- 
spins is possible 0, whose properties we now review. 
Once again, we concentrate on the case where the bosons 
are strongly interacting ([/bb = 8) and all the cases where 
A^f -|-iVb = L. First, for small Uhf, the full system is com- 
posed of two weakly interacting Luttinger liquids, one for 
each species. Green functions decay algebraically (Fig.jSj 
L^bf < 2) and the winding numbers Wh and W{ are not 
correlated (Fig. [SI left). 

For large enough Uhf the system enters a PSDW phase: 
The interactions forbid double occupancy of a site by 
two bosons or by a fermion-boson pair and a pseudo- 
spin along the z-axis on site i may be defined as erf = 
n,^ — = ±1. The model can then be mapped on the 
XXZ Hamiltonian 



Ti-xxz = Jx 
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where J^y ^ -t^ /Um and t^il/UM - l/J/bb). The 
magnetic correlations along the xy axes correspond to 
particle- hole correlations Ghf{x) as af oc b\fi. Bosonic 
and fermionic Green functions decay exponentially as the 
hopping of individual particles is energetically suppressed 
(Fig. [5l Uhf > 3). The hopping of bosons and fermions 
and, consequently, the winding numbers are strongly 
anti-correlated Wh = —Wf (Fig. [HI right) which leads 
to pb = Pf = Pbf/4 and p'^^ = 0. 
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FIG. 5: Green functions for several Um- At low interaction, 
there are superfluid bosons and mostly free fermions (Lut- 
tinger liquids). At large interaction Uhf > 3, the system enters 
the PSDW regime, individual Gb and Gf decay exponentially 
while Gbf remains quasi long ranged. 
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FIG. 6: Histograms of the winding numbers J\f{Wh,Wf). 
Same parameters as in Fig. [5] In a low interaction Luttinger 
liquid phase {Uhf = 1) Wh and W{ are not correlated and A/" is 
isotropic (left panel). With increasing Uhf the winding num- 
bers become correlated (middle panel), while in the PSDW 
phase {Uhf = 4), A/" is strongly anisotropic and only allows 
the cases where Wh ~ —Wt (right panel). 



Along the z axis, one defines a spin correlation func- 
tion Cpsdw(a;) = (cf crf_|_^). Depending on the value of Jz, 
the system favors antiferromagnetic (AF), for Uhf < Uhh, 
or ferromagnetic (FM), for Uhf > Uhh, correlations along 
the z-axis. For Uhf > 2Uhh, \Jz\ > \Jxy\, the dominant 
correlations are ferromagnetic along the z-axis which 
leads to phase separation for the bosons and fermions. 
In all other cases, the xy term is dominant, i.e., Gbf is 
the correlation function that has the slowest decay. 

We concentrate on this latter case. In the AF regime 
along the z direction, we observe the characteristic oscil- 
lations around zero, the FM correlations decreasing much 
more rapidly than the AF ones (see Fig. [7]). The pres- 
ence of defects in the AF due to the imbalance between 
the populations of up and down spins leads to oscillations 
with beating where Gpsdw(a;) behaves like cos{2TTNfx/L) 
multiplied by a power law envelope. Density correla- 
tion functions exhibit similar beating oscillations, rem- 
iniscent of the solitonic excitations observed in the ex- 
tended Hubbard model [l^. The commensurate oscilla- 
tions in cos(7ra;) are recovered only for iVf = iVb = L/2 
(see Fig. [T]). 
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FIG. 7: Pseudo spin correlation function Cpsdw(2;) in the case 
A*'f — Nh — L/2 for the AF (circles), pure-a;y (triangles) and 
FM (diamonds) regimes. Inset: Cpsdw(a;) in the CSF phase for 
different values of il7bf- The AF signal decreases as expected. 



On the contrary, in the FM regime, we see no sign of 
oscillations and observe only a slow FM relaxation. For 
the special case where Jz = {Um — Uhh) the system is 
mapped on a pure xy pseudo-spin Hamiltonian. The AF 
and FM correlation functions are expected to decay with 
the same power law [l6| and we observe the combined 
effect of both types of correlation (see Fig. [7]). 

In the SFb phase we have already seen that xy corre- 
lations are present (Fig.|4]). It is more difficult to observe 
the PSDW correlations along the z axis. To obtain the 
SFb phase, one needs Uhi > Uhh for the Mott phase to 
be destroyed. This means that the filled layer is in the 
FM regime. However, it is difficult to distinguish the FM 
signal due to the filled layer from the signal due to the 
superfiuid boson density-density correlations which are 
expected to show some kind of short range relaxation. 
On the contrary, an AF signal would be due to the filled 
layer only, but should be small because we are in the re- 
gion where the dominant effects are FM. We observe such 
a weak AF signal with a decreasing amplitude when Uhf 
increases, as expected (Fig. [3 inset). 

Finally the interpretation of the SFb as composed of 
superfiuid bosons on a filled layer in a pseudo-spin phase 
is partly supported by our data. We observe, as expected, 
that the leading correlations are the bosonic Green func- 
tions followed by the composite Green function. We also 
found that the measure of the pseudo-spins correlations 
along z are consistent with the PSDW phase. However, 
the decay of the fermionic Green function is algebraic 
which is consistent with the non-dominant PSDW corre- 
lations. The present algorithm gives access to one- and 
two-body correlation functions and topological informa- 
tion via the winding numbers. It does not give informa- 
tion about other A^— point correlation functions or what 
are the dominant correlations. 

In conclusion, we have studied the effect of adding 
fermions on a boson Mott insulator and varying the in- 



teraction strength between fermions and bosons. We ob- 
served that, when the Mott insulator is destroyed, the 
system enters a superfiuid phase where the movement of 
bosons and fermions is correlated. The similarities be- 
tween this phase and a previously observed pseudo-spin 
phase led us to describe it, approximately, as a phase 
displaying simultaneously traditional superfluidity and 
PSDW. 

Even in non symmetric cases, where the number of 
bosons and fermions is not equal, it is possible to observe 
some of the characteristics of the original pseudo-spin 
phase as long as interactions are large. These results 
lead us to conclude that anticorrelated motion of bosons 
and fermions is a robust property and appears generally 
in strongly interacting mixtures of bosons and fermions. 

The authors thank U. SchoUwock, B. Svistunov, and 
M. Troyer for useful discussions. We acknowledge finan- 
cial support from a grant from the CNRS (France) PICS 
18796 and from the Swiss National Science Foundation. 
Part of the calculations were carried out on the Hreidar- 
cluster at ETH Zurich. 



[1] 
[2] 



[3] 



[4] 
[5] 

[6] 
[7] 
[8] 

[9] 
[10] 

[11] 
[12] 



M. Greiner, O. Mandel, T. Esslinger, T.W. Hansch, and 
I. Bloch Nature 415 39 (2002). 

J.E. Lye, L. Fallani, M. Modugno, D.S. Wiersma, C. Fort, 
and M. Inguscio Phys. Rev. Lett. 95, 070401 (2005). 
D. Clement, A.F. Varon, M. Hugbart, J.A. Ret- 
ter, P. Bouyer, L. Sanchez-Palencia, D.M. Gangardt, 

G. V. Shlyapnikov, and A. Aspect, Phys. Rev. Lett. 95, 
170409 (2005). 

F. Schreck, L. Khaykovich, K. L. Corwin, G. Fer- 
rari, T. Bourdel, J. CubizoUes, and C. Salomon, Phys. 
Rev. Lett. 87, 080403 (2001). Z. Hadzibabic, S. Gupta, 
C.A. Stanm C.H. Schunck, M.W. Zwierlein, K. Dieck- 
mann, and W. Ketterle, Phys. Rev. Lett. 91, 160401 
(2003). 

M.W. Zwierlein, C.H. Schunck, A. Schirotzek, W. Ket- 
terle, Nature 442, 54 (2006). 

K.Giinter, T. Stoferle, H. Moritz, M. Kohl, and 
T. Esslinger, Phys. Rev. Lett. 96, 180402 (2006). S. Os- 
pelkaus, C. Ospelkaus, O. Wille, M. Succo, P. Ernst, 
K. Sengstock, and K. Bongs Phys. Rev. Lett. 96, 
180403 (2006). 

A.B. Kuklov and B.V. Svistunov, Phys. Rev. Lett. 90, 
100401 (2003). 

M. Lewenstein, L. Santos, M.A. Baranov, and 

H. Fehrmann, Phys. Rev. Lett. 92, 050401 (2004). 

L. Mathey, D.W. Wang, W. Hofstetter, M.D. Lukin, 
and E. Demler, Phys. Rev. Lett. 93, 120404 (2004); 
L. Mathey and D.W. Wang, Phys. Rev. A 75, 013612 
(2007). 

P. Sengupta and L.P. Pryadko"cond-mat /0512241 1 (2005) . 
L. Pollet, M. Troyer, K. Van Houcke, and S.M.A. Rom- 
bouts, Phys. Rev. Lett. 96, 190402 (2006). 
L. Pollet, C . KoUath, U. SchoUwock, and M. Troyer, 
|cond-mat/06 09604 021602(R) (2006). 
S.M.A. Rombouts, K. Van Houcke, and L. Pollet, 



5 



Phys. Rev. Lett. 96, 180603 (2006). K. Van Houcke, 
S.M.A. Rombouts, and L. PoUet, Phys. Rev. E 73, [15] 
056703 (2006). 

[13] L. Pollet, Ph. D. Thesis, Universiteit Gent (2005 ). See [16] 
|http : //www.nustruc .ugent .be/thesislode .pdf I 

[14] These two regions are essentially the same as the two 



Luttinger's liquid regions (LLa and LLb) found in [T^ 
F. Hebert, G.G. Batrouni, and R.T. Scalettar, Phys. Rev. 
A 71, 063609 (2005), [cond-mat/0412037| 
T. Giamarchi, Quantum Physics in One Dimension, Ox- 
ford (2004). 



